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ABSTRACT 

We present Cosmological models with modified Chaplygin gas (MCG) in the 
framework of Horava-Lifshitz (HL) theory of gravity both with and without detailed 
balance. The equation of state (EOS) for a MCG contains three unknown parameters 
namely A, a, B. The allowed values of some of these parameters of the EOS are 
determined using the recent astrophysical and cosmological observational data. Using 
observational data from H(z) — z, BAO peak parameter, CMB shift parameter we 
study cosmologies in detailed-balance and beyond detailed-balance scenario. In this 
paper we take up the beyond detailed-balance scenario in totality and contribution of 
dark radiation in the case of detailed-balance scenario on the parameters of the EOS. 
We explore the effect of dark radiation on the whole range the of effective neutrino 
parameter to constrain matter contributing parameter B in both the detailed-balance 
and the beyond-detailed balance scenario. It has been observed that greater the dark 
radiation less the matter contribution in the MCG in both the scenario considered 
here. In order to check the validity of beyond detailed balance scenario we plot su- 
pernovae magnitudes (/z) with redshift of Union2 data and then the variation of state 
parameter with redshift is studied. It has been observed that beyond detailed balance 
scenario is equally suitable in HL gravity with MCG. 
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1 INTRODUCTION 

The big-bang cosmology has become the standard model for 
cosmology which accommodates a beginning of the Universe 
at some finite past. The discovery of CMBR (Penzias, Wil- 
son 1965; Dicke et al. 1965) supports such model of the uni- 
verse. However, big-bang cosmology based on perfect fluid 
assumption fails to account some of the observed facts both 
in the early and at late universe. The standard Big-bang 
model is known to have several limitations; for instance, (i) 
Horizon problem (ii) Flatness problem (iii) the singularity 
problem, to name a few. It is known that the above prob- 
lems can be resolved by invoking a phase of inflation at a 
very early epoch. Most of these problems have, however, 
been resolved by invoking inflation (Guth 1981; Linde 1982; 
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Albrecht & Steinhardt 1982; Sato 1981) in the semiclassi- 
cal theory of gravity. On the other hand recent observa- 
tions predict that our universe is passing through a phase of 
acceleration (Riess et al. 1998). This phase of acceleration 
is believed to be a late time phase of the universe and it 
comes out that such a phase cannot be accommodated in 
the general theory of relativity with the usual matter fields 
in the standard model of particle physics. Since a universe 
with inflation should give satisfactory explanation of what 
happens close to the Planck era, it is also necessary to con- 
sider a satisfactory theory which is valid near that epoch. 
It may be pointed out here that a quantum gravity effect 
becomes important at the Planck time. But a consistent 
theory of quantum gravity is yet to emerge. In this direction 
superstring theory may be considered as one of the promis- 
ing candidate of quantum theory of gravity. Cosmological 
models are also proposed in Loop Quantum Gravity (LQG) 
(Bojowald 2001) which avoids initial singularity. However, a 
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proper description of time evolution of quantum space-time 
in the LQG is not satisfactory. Several attempts have been 
made in the recent past to achieve a complete quantum grav- 
itational theory (UV complete theory). Among many such 
attempts, Horava-Lifshitz (henceforth, HL) theory of grav- 
ity appears to be interesting. The success of the Lifshitz 
theory in solid state physics motivated Horava to propose a 
theory of gravity, often called Horava-Lifshitz (HL) gravity 
(Horava 2009) which may be important to explore a vi- 
able cosmological model. In the ultraviolet (UV) limit, HL 
gravity has a Lifshitz-like anisotropic scaling as t — > Ft and 
x % — > lx l , between space and time, characterized by the dy- 
namical critical exponent z — 3 and thus breaks the Lorentz 
invariance; while in the infra-red (IR) limit, the scale reduces 
to z = 1. So, it is expected that it may reduce to classical 
general relativistic theory of gravity in the low energy limit. 
The Friedmann equation gets modified by a -\ term (Lu, 
Mei & Pope 2009; Calcagni 2009; Kiritsis & Kofinas 2009), 
where a is the scale factor in a non-flat universe in the HL- 
gravity. 

In the original HL gravity, Horava assumed two con- 
ditions: detailed balance and projectibility. More recently, 
Sotiriou, Visser and Weifurtner (SVW) (Sotiriou, Visser & 
Weinfurtner 2009) , proposed a general HL theory with pro- 
jectability but without detailed-balance conditions. For a 
spatially curved Friedmann- Robertson- Walker universe, the 
SVW generalization yields an extra -term that modifies 
the coefficient of the -\ term in the Friedmann equation as 
compared to the HL theory. Therefore, it is important to 
look for cosmological models in Horava gravity considering 
projectibility with and without detailed-balance. 

In the HL-gravity, the initial bigbang singularity may 
not arise due to the presence of higher order terms in the 
spatial curvatures Rij [4]. There are many such novel fea- 
tures of HL gravity for which it is worth to explore different 
aspects of observed universe. A volume of literature in the 
framework of the HL gravity appeared containing the study 
of gravitational wave production (Mukohyama et al. 2009; 
Park 2009; Myung 2009), perturbation spectrum (Gao et. al. 
2009; Cai and Zhang 2009; Wang and Maartens 2010), black- 
hole properties (Danielsson and Thorlacius 2009; Cai, Cao 
& Ohta 2009; Kehagias & Sfetsos 2009), dark energy phe- 
nomenology (Park 2010; Chaichian et. al. 2010), the prob- 
lems of determining observational constraints in the theory 
(Dutta & Saridakis 2010), astrophysical phenomenology 
(Kim et. al. 2009; Harko 2009; Iorio and Ruggiero 2009), 
thermodynamical properties etc. (Wang & Wu 2009; Cai, 
Cao & Ohta 2009). Though there exists foundational and 
conceptual issues of Horava-Lifshitz gravity and its associ- 
ated cosmology, cosmological scenario have been examined 
with generalised Chaplygin gas (GCG) (Ali et. al. 2010). 
GCG being an exotic matter may be useful to address the 
recent acceleration of the universe. One of the characteristic 
features of the GCG is that it behaves as a pressureless fluid 
at the early stage of the evolution of the universe, and at a 
later stage it behaves as a cosmological constant. Recently, 
a modified form of Chaplygin gas is also considered exten- 
sively in cosmology (Liu & Lu 2005; Thakur, Ghose & Paul 
2009). The modified Chaplygin gas (MCG) is more general 
and contains three free parameters. The idea is to interpolate 
states of standard fluids at high pressures and at high en- 
ergy densities to a constant negative pressure at low energy 



densities (Debnath, Banerjee and Chakraborty 2004) . In the 
present work we explore cosmological models with MCG in 
the framework of HL gravity and determine the range of pa- 
rameters of MCG from recent cosmological observations. We 
examine the effect of effective neutrino parameter on both 
detailed-balance and beyond-detailed balance scenario. Here 
compatibility of beyond detailed balance scenario is also ex- 
plored in details using recent observational data in HL grav- 
ity with MCG . The objective of the paper is to determine 
the limits of the unknown EOS parameters using the obser- 
vational data. The equation of state parameter of the total 
cosmic fluid defined by w(z) = will be evaluated and 

examined at different values of redshift parameters. Com- 
paring the supernovae magnitudes (/i) vs. redshift (z) with 
Union2 data we test the viability of beyond-detailed balance 
scenario. The suitability of the model is also examined using 
w(z) vs. z plot in our model. 

The paper is organized as follows: In sec. 2, we present 
the basic equations for Horava-Lifshitz cosmology and ob- 
tain the Friedmann equations for detailed balance and be- 
yond detailed balance conditions. In sec. 3, the energy den- 
sity and EOS for MCG is presented. In sec. 4 , the con- 
straints on detailed-balance condition and beyond detailed 
balance condition from the observations is presented. In sec. 
5, numerical analysis to determine constraints on EOS pa- 
rameters are obtained for detailed balance. In sec. 6, numer- 
ical analysis to determine constraints on EOS parameters 
are obtained for beyond detailed balance scenario. In sec. 7, 
the viability of MCG in HL gravity is discussed. Finally, in 
sec. 8, we summarize the result. 



2 HORAVA-LIFSHITZ COSMOLOGY 

In Horava-Lifshitz gravity (Calcagni 2009; Kiritsis & Kofinas 
2009), it is convenient to use the four-dimensional space-time 
metric of the Arnowitt-Deser-Misner (ADM) decomposition 
form which is given by 



ds 2 = -N 2 dt 2 + giJ (dx* + N l dt)(dx j + N j dt) 



(1) 



where the basic variables are lapse function N, shift vector 
Ni, and the spatial metric gtj. The scaling transformation 
of the co-ordinates reads: t — > l 3 t and x x —* lx l . The shift 
N l and the 3d spatial metric gij depend both on the time 
coordinate t and the spatial coordinate x 1 , the lapse N is as- 
sumed to depend on time only. This condition on the lapse is 
called the projectibility condition. The action of HL gravity 
consists of kinetic and potential terms as follows: 



(2) 



S g = S K + Sv = J dt d'x^g N(L K + L v ). 
The kinetic term is given by 



-/ 



S K = I dtd 3 x^N 



2(K ij K ij - XK 2 



(3) 



where Ki 



^ — V '^- VjJV '^ is the extrinsic curvature and 



dot represents derivative with respect to time (£). 



2.1 Detailed Balance Condition and Projectibility 

The symmetry property of the Lagrangian Lv , employed in 
the gravitational action drastically reduces the number of in- 
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variants which one should actually consider in the action to 
begin with (Horava 2009). The above symmetry is known as 
detailed balance which follows from condensed matter sys- 
tems and requires that the Lagrangian Lv should be deriv- 
able from a superpotential W (Kiritsis & Kofinas 2009). 
Under the detailed balance condition the total action of HL 
gravity is given by 



dtd 3 x^/g~N 



■2il\,;K ■ - \K 2 ) k 2 C 1] C 13 

K 2 + 2iJ 4 

K 2 ne ijk RiiVjR l k K 2 p?RijR i3 ' 



1 2 2 

k p, 



8(3A- 1) 



(l-4A)i? 2 



+ AR- 3A 2 



where 



ijk 



C 13 = -V t 
y/9 



(4) 



(5) 



is known as the Cotton tensor, and the covariant derivatives 
are determined with respect to the spatial metric (gij), e t3k 
is a totally antisymmetric unit tensor, A is a dimensionless 
constant and the variables k, u> and p are constants . 

In the above gravitational action to include matter com- 
ponents one needs to add a cosmological stress-energy ten- 
sor to the gravitational field equations, that recovers the 
usual general relativity formulation in the low-energy limit 
(Sotiriou, Visser & Weinfurtner 2009; Chaichian et. al. 2010; 
Carloni, Elizalde & Silva 2009). The matter-tensor is a hy- 
drodynamical approximation that leads to the existence of 
energy density (p m ) and pressure (p m ) in the Friedmann 
equation, where p m represents the total matter energy den- 
sity, that accounts for both the baryonic pb as well as that 
of the dark matter p dm , including the normal matter (where 
p m represents pressure). 

Horava obtained the gravitational action assuming that 
the lapse function is just a function of time i.e., N = N(t). 



Here we use FRW metric with N 
with 



dr 2 



1 - Kr 2 



+ r 2 dD.2 



I, gij =a 2 (t)^ij,N l =0 



(6) 



where K = — 1, 1, 0, corresponds to open, close and flat uni- 
verse respectively. By varying N and gij in the gravitational 
action (4), one obtains the following field equations: 

H 2 = TTTTTT" —(Pm+Pr) 



6(3A - 1) 

K 2 

+ 6(3A- 1 

n 4 p 2 AK 
8(3A- l) 2 a 2 



3k 2 p 2 K 2 3n 2 p 2 A 2 
8(3A - l)a 4 + 8(3A - 1) 



(7) 



H + 



3H 2 



4(3A - 1) 
k 2 

4(3A-1) 

K 4 p 2 AK 
"16(3A- l) 2 a 2 



{PmiOm, + Pr^r) 



_ 3*VA 2 
8(3A - l)a 4 8(3A - 1) 



K 2 p 2 K 2 



(8) 



where H — f- In the above field equations the term 
proportional to a~ 4 may be considered as the usual 
" dark radiation term" , characteristics of the HL Cosmology 
(Calcagni 2009; Kiritsis & Kofinas 2009) and the constant 
term is identified with the usual cosmological constant. The 
conservation equation for matter is: 



pm + 3H(p m + p m ) = 0, 

and that of radiation is: 

Pr + 3H(p r +p r ) = 0, 

where we denote 



K 



8(3A - 1; 



16tt(3A- 1) : 
5T = 1. 



2.2 



32tt' 

Beyond Detailed Balance Condition with 
Projectibility 



(9) 
(10) 

(11) 
(12) 
(13) 



As it is not known with certainty whether the detailed bal- 
ance condition is enough for extracting whole information 
of Horava-Lifshitz gravity (Calcagni 2009; Kiritsis & Kofi- 
nas 2009) or it is necessary to do something with this balance 
condition, one can investigate cosmological scenario in the 
HL gravity relaxing the detailed balance condition. In this 
subsection we discuss the cosmology of HL gravity in the 
pressence of modified Chaplygin gas (MCG), baryon, radia- 
tion and dark radiation without detailed-balance. The aim 
of the paper is to look for the effects of the dark-radiation 
on the parameters of the MCG model. The Friedmann equa- 
tions in this case can be written as (Sotiriou, Visser & We- 
infurtner 2009; Carloni, Elizalde & Silva 2009; Bogdanos & 
Sarikdakis 2010; Charmousis et. al. 2009; Leon & Saridakis 
2009): 

2<7 



H z = 



(3A-1) 



{pm + Pr) 



+ 



(3A-1) 
3(3A - l)a 2 



£l 

6 



+ 



a 3 K 2 a 4 K 



6a 4 



6a 6 



(14) 



H+ 2 H - (3A - T) 



+- 



(3A - 1) 
0-2K 



{pm^m + prOJr 



0\ U3K O^K 

"~6~ + 18a 4 + "6^6" 



6(3A- l)a 2 ' 

where ao = k 2 /12. 

The dimensionless parameters are given by: 

_ 60-0 
o cosmo - g ^ (3A _ x) , 

ct 2 = -3(3A- 1), 



G a 



6(70 

167r' 



(15) 

(16) 
(17) 
(18) 
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where 01 < and 04 > 0. In the case of the detailed balance 
case, in the IR limit (A = 1), the two parameters G CO smo 
and Ggrav coincides. 



3 EOS FOR MODIFIED CHAPLYGIN GAS 

The equation of state for Generalized Chaplygin Gas (GCG) 
(Billic, Tupper & Viollic 2001; Bento, Bertolami & Sen 
2002) is given by 



A_ 



(19) 



with < q < 1. In the above original Chaplygin gas corre- 
sponds to a = 1 (Chaplygin 1904). It may be pointed out 
here that Chaplygin introduced his equation of state (Chap- 
lygin 1904) to study the lifting force on a plane wing in 
aerodynamics. Chaplygin's equation of state (19) has raised 
recently a renewed interest (Bazeia 1999; Jackiw & Poly- 
chronakos 1999) because of its many remarkable and, in 
some sense, intriguingly unique features. It is interesting to 
note that it has amazing connection with string theory: it 
can be obtained from the Nambu-Goto action for d-branes 
moving in a (d+2)-dimensional spacetime in the light-cone 
parametrization (Bordemann & Hoppe 1993). However, in 
cosmological context generalized form of Chaplygin gas may 
be useful to describe the observed universe. GCG has two 
free parameters A (positive), a. Recently a further mod- 
ification of GCG has been proposed in the framework of 
cosmology (Liu & Lu 2005). The modified Chaplygin gas 
(MCG) is more general and it contains one more free pa- 
rameter (B). The model is consistent with (i) Gravitational 
lensing test (Silva & Bertolami 2003; Dev & Alcaniz 2004) 
and (ii) Gamma-ray bursts (Bertolami & Silva 2006). The 
equation of state for the MCG is given by: 



p = Bp- 



A 



(20) 



where A, B, a are arbitrary constants to be determined from 
observation for model building with < a < 1. 

The energy conservation equation for the MCG is: 



p c + 3H(p c +Pc) =0 



(21) 



where p c and p c correspond to energy density and pressure 
of MCG respectively. Using eq. (20) in eq.(21) we obtain: 



Pc = 



A C 
1 + B 



3n 



(22) 



where C is an arbitrary constant and we denote (1 + B)(l + 
a) — n. Equation (22) can be rewritten as 



pc = po 



A s + 



As 



(23) 



where we denote As = Jt+i with = , z is redshift 



l + z' 



parameter and we choose a = 1 for convenience. MCG 
reduces to GCG model when we set B — in the above 
equation. 



4 OBSERVATIONAL CONSTRAINTS ON EOS 
PARAMETERS 

In general theory of relativity cosmological models with 
MCG has been studied and the constraints on EOS param- 
eters for viable cosmologies are determined using observa- 
tional data (Thakur, Ghose & Paul 2009; Jianbo et al. 2010). 
The EOS parameters of MCG will be explored here for vi- 
able cosmologies in the framework of Horava-Lifshitz gravity 
using the recent observational data. For this we have taken 
up data from Observed Hubble Data (OHD), BAO peak 
parameter and CMB shift parameter. 



4.1 Constraints Obtained from Detailed Balance 

In this case, using eqs. (7) and (8), the Friedmann equations 
can be rewritten as: 



= "I" + * + 2A^ + 2 



K_ 

^2' 



tt 3 rr'2 A 1 x ( K 2 3A 

H+ 2 H = " 4 - G ^ + 3^ ) -( v ^ I "~ 
K 



(24) 



(25) 



Let us define the following dimensionless density pa- 
rameters: 

(i) for matter component 



8ttG 
1 ~ 3/P Pl 



(ii) for curvature 



K 



H 2 a 2 

(iii) for cosmological constant 



n o ee 



2HI 



(26) 



(27) 



(28) 



We define another dimensionless parameter for expan- 
sion rate as: 



E(z) 



H(z) 
H ' 



(29) 



Using the above definition of parameters, the Friedmann 
equation now can be rewritten as: 

e 2 {z) = n bo (i + z) 3 + n C oF(z) + n ro (i + z) 4 



+n Ko {i + z) 2 + o + 



^Ko(l + 2) 



4ft 



where 



F(z) 



A s + 



I- As 



3(l+B)(l+a) 



(30) 



(31) 



Let us assume E{z = 0) = 1 at the present epoch, which 
leads to 



Viho + Qco + + Q.K0 + flo + 



4fi D 



(32) 
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where Qbo, ^ co , S7ro, fi/fo represent the present day baryon, 
MCG, radiation and curvature energy density parameters 
respectively. Here fl is the energy density associated with 
the cosmological constant. The last term in eq. (32) corre- 
sponds to dark radiation, which is a characteristic feature 
of the Horava-Lifshitz theory of gravity. The dark radiation 
component may be important during nucleosynthesis. Thus 
a suitable bound from Big Bang Nucleosynthesis (hence- 
forth, BBN) may be incorporated in the above EOS. Using 
the upper limit on the total amount of Horava-Lifshitz dark 
radiation that is permitted during BBN era is expressed by 
the parameter AN V which represents the effective neutrino 
species (Hagiwara et. al. 2002; Olive, Steigman & Walker 
2000). We obtain the following constraint equation (Dutta 
& Saridakis 2010): 



o 2 

= 0.135AJV„n ro 

4f2 



(33) 



The BBN upper limit on AiV„ is -1.7 AiV„ ^ 2.0, is taken 
from the Refs. (Olive, Steigman & Walker 2000; Steigman 
2006). A negative value of AN V is usually associated with 
models involving decay of a massive particles which we do 
not consider here. Again AN V — 0, which corresponds to the 
zero curvature scenario will be excluded also. It is because of 
the fact that the Horava-Lifshitz cosmology with zero curva- 
ture becomes indistinguishable from ACDM. The curvature 
in dynamical dark energy models are important, neglect- 
ing the curvature term impose a serious problem (Clark- 
son, Cortes & Bassett 2007; Virey et. al. 2008). Therefore, 
we consider the limiting values for AN V which satisfies the 
bound < AA^ sC 2.0. 

The numerical analysis taken up here contains nine pa- 
rameters, these are namely, Qbo, f^co, f!r», f2 , AN V , 
H , As, B, a. As the number of unknowns are more than 
the number of equations we fix some of the parameters us- 
ing the best-fit values from 7 year WMAP data (Komatsu 
et. al. 2010). The fixed parameters are Q mo (= Qto + fico), 
flbo, H , flro and the corresponding values of the param- 
eters are chosen as follows : Q. mo = 0.27, Qbo = 0.04, 
Ho = 71AKm/sec/Mpc, Q ro = 8.14 * 10" 5 . Therefore, one 
can have now only six free parameters to be determined 
which are Qko, fl , As, B, a, AN„. Using eq. (33) in eq. 
(32) one obtains 

Q (K, AN U , As, a) = 1 - n mo - (1 - 0.135AJV I ,)a. o 



-0.73(K)y/AN^y/Q ro - n mo « ro - fip (34) 
Q.Ko{AN v , As, a) = y/0.bAAN„n ro n o {K, AN v ,A s ,a) (35) 

which may be employed for a close or in an open universe 
depending on the values of K. Now, it reduces to four free 
parameters, namely, As, B, a, AN V . To determine the effect 
of dark radiation on the constraints of the parameters of 
the MCG in detailed-balance scenario, we took two extreme 
values of a (a=0.999, 0.001) satisfying < a < 1 for two 
extreme values of AN V (0.01, 2.0) in both close and open 
universe respectively. In this case each of these values of 
q and AJV„ determined the best-fit values of the rest two 
parameters (i.e., As, B). Thereafter at the extreme values 
of AN V for two extreme values of a we plot contours for 
the parameters As, B at different confidence levels. From 
the contours of As, B drawn at different values of a and 
AJV„ we determine the permissible range of values of the 



B-parameter for the MCG in HL gravity in the framework 
of open or closed universe. 

4.2 Constraints Obtained from Beyond-Detailed 
Balance 

In beyond-detailed balance scenario using eqs. (14)-(15), the 
Friedmann's equations can be rewritten as: 

8nG, 



H 



-(Pb + Pc + Pr 



6 + 6a 4 



K_ 
775" 



(36) 



h + Ih> 



-4nG(p c + -p r ) 



3 
2 
K 
2^ 



' 6 18a 4 6a 6 



(37) 



Now one can re-write the dimensionless Hubble param- 
eter as follows: 



E 2 {z) = Q bo (l + z) 3 +n co F(z) + n ro (l + z) 4 + tt Ko (l + z) 2 



where 



F{z) = 



A s + 



+{n 1 + n s (i + z) 4 + n 4 (i + zf} 

l-As 



(38) 



(39) 



a 3(l+S)(l+a) 

The dimensionless parameters namely, fii, SI3, Q4 are re- 
lated to the model parameters G\ , 03 , 04 as follows: 

2r>2 



n 3 = 



Q4 



aiH Vt Ko 
6 

o/SIko 



(40) 

(41) 
(42) 



At the present epoch E(z = 0) = 1, which leads to 

fifm + !2 co + H r „ + Q.Ko + fii + ^3 + «4 = 1 (43) 

In the above equations Q4 is required to be a posi- 
tive quantity in order that the Hubble parameter remains 
positive at all redshifts, also the gravitational perturbations 
(Sotiriou, Visser & Weinfurtner 2009; Bogdanos & Sarik- 
dakis 2010) demand the same. For the conveneince of our 
analysis ^3 is assumed positive definite. 

The above equation is a constraint equation for this 
analysis and we use it to replace f2i in terms of other pa- 
rameters in our analysis. Following the procedure adopted 
in Ref. (Dutta & Saridakis 2010) for AN V , we consider 
the upper limit of dark radiation beyond standard model 
at the BBN. Consequently, the following constraints at the 
time of BBN emerged (z = zbbn) (Hagiwara et. al. 2002; 
Olive, Steigman & Walker 2000; Steigman 2006; Malaney 
and Mathews 1993) : 

3 + 4 (1 + z 2 BBN f = Q, 3max = 0.135AA^n r o. (44) 
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where the 03 represent the usual dark radiation and rep- 
resents a kinetic-like component (a quintessence field dom- 
inated by kinetic energy) (Joyce 1997; Joyce and Prokopec 
1998). The above equation will be used to replace fi4 in 
terms of other parameters in the analysis. For simplicity we 
define 

f* = (45) 

where Q^max is the upper limit on £1$. This will help us to 
express ^3 in terms of other parameters. 

Following the detailed-balance scenario we consider 
AN V here to satisfy the bound < AN V < 2.0, following the 
importance of curvature in dark energy models and treating 
Qko as a free parameter (Clarkson, Cortes & Bassett 2007; 
Virey et. al. 2008). 

To sum up, in the numerical analysis taken up here, 
the following parameters, Qbo, ^co, llro, Qko, fii, ^3, O4, 
AN V , H , As, B, a, /3 are involved. We fix some of the 
parameters using the best-fit values from 7 year WMAP data 
(Komatsu et. al. 2010). The fixed parameters are Q. mo {= 
Q,bo + Sim), fi&o, H , fl ro and the corresponding values of 
the parameters are chosen as follows : fl mo = 0.27, flbo = 
0.04, H = 71AKm/sec/Mpc, fi ro = 8.14 * 10 -5 . Using 
the constraint eqs. (36) - (42) one can replace Qi, O3 , f24 
in terms of the other six free parameters for the numerical 
analysis. Therefore one can have only six free parameters to 
be determined which are Qko, As, B, a, j3, AN V . 

To determine the constraints on the parameters of the 
MCG in beyond detailed-balace scenario, we consider three 
values of a satisfying < a < 1 (a=0.999, 0.500, 0.001) 
and determine the best-fit values for the rest five parame- 
ters (i.e., As, B, (5, Qko, AN v ). Thereafter, at the best-fit 
values of AN„, f3, Qko for those three values of a we plot 
2d contours for the parameters As, B at different confidence 
levels. The contours of As, B drawn at different values of 
a in turn determine the permissible range of values of the 
B-parameter for the MCG in HL gravity in the framework 
of beyond detailed balance scenario. 

To examine the effect of dark radiation (i.e., effective 
neutrino parameter) on the constraints on the parameters 
of the MCG, we took two extreme values of a (a=0.999, 
0.001) satisfying ^ a ^ 1 for two extreme values of AN V 
(0.01, 2.0). In this case each of these values of a , AN V 
determines the best-fit values of the rest four parameters 
(i.e., As, B, (3, Qko)- Thereafter, at the extreme values of 
AJV„ for two extreme values of a we plot 2d contours for 
the parameters As, B for the best-fitted values of /3, O.K0 
at different confidence levels. From the contours of As, B 
drawn at different values of a and AN V we determine the 
permissible range of values of the B-parameter for the MCG 
in HL gravity in the framework of beyond detailed balance 
scenario. We note that the range of values of B is narrower 
due to the effect of effective neutrino parameter on B. 



set for (H-z) data has been used along with BAO peak pa- 
rameter and CMB shift parameter. Chi-square minimisation 
technique has been used here to determine the limiting val- 
ues of the EOS parameters in the next subsections. 



5.1 (H-z) data as a tool for constraining 

The best-fitted parameters of the model considered here can 
be obtained by minimising the entity chi-square which is 
defined as 



XOHD(H ,A s ,B,a,AN v ,z) = 
{H{Ho,A s ,B,a,AN v ,z) - H bs{z)) 2 



(46) 
(47) 



where H b B (z) is the observed Hubble parameter at red- 
shift(z) and a\ is the associated error with that particular 
observation. The Hubble parameter is given by 



H(z) = H E(z) 
where 



E(z) 



{Q bo (l + zf + n co F{ Z ) + Q ro {l + . 



+ V K o(1+z) 2 +[{}o+ Qk ° {1 + Z) 



and denoting 



F{z) = 



A s + 



As 



a 3(l+S)(l+a) 



(48) 



) 1/2 (49) 



(50) 



Here H(z) — z data is taken from Stern Data analysis (Stern 
et. al. 2010). There are 12 data points of H(z) at redshift z 
used to constrain the MCG model. 



5.2 BAO peak parameter as a tool for 
constraining 

A model independent BAO (Baryon Acoustic Oscillation) 
peak parameter can be defined for low redshift (zi) mea- 
surements as: 



A- 



[E(zi)]t 



T z l dz 
JO E(z) 



2 
3 



Zl 



(51) 



where O m is the matter density parameter for the Universe. 
For a detailed description of the above defined parameter 
and related approximations reader is referred to ( Eisenstein 
et. al. 2005). The chi square function can be defined as: 



2 _ {A- 0.469) 1 

Xbao — 



(0.017)' 



(52) 



where we have used the measured value for A (0.469±. 0.017) 
as obtained by ( Eisenstein et. al. 2005) from the SDSS data 
for LRG (Luminous Red Galaxies) survey. 



5 NUMERICAL ANALYSIS TO DETERMINE 
CONSTRAINTS ON THE EOS 
PARAMETERS IN DETAILED BALANCE 
SCENARIO 

In this section we use three sets of different observational 
data to constrain the parameters of the MCG. Stern data 



5.3 CMB Shift Parameter as a tool for 
constraining 

Here the CMB shift parameter is defined as 

rz, ° dz 



f 

Jo 



E(z) 



(53) 
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Figure 1. Constraints for closed universe for a = 0.999 from 
OHD+SDSS+CMB Shift data 68.3%(Solid) 95.4% (Dashed) and 
99.73% (Dotted) contours. 



a = .001 
AN^ 2.0 




Figure 2. Constraints for closed universe for a = 0.001 from 
OHD+SDSS+CMB Shift data 68.3%(Solid) 95.4% (Dashed) and 
99.73% (Dotted) contours. 



Model 


B 


A s 


AN V 


a = 0.999 


0.003745 


0.062817 


0.232994 


a = 0.500 


0.016592 


0.110548 


0.099996 


a = 0.001 


0.006192 


0.052076 


0.807051 



Table 1. Best-fit values for MCG: K = 1 

where zi B is the z at last scattering. The WMAP7 data gives 
us R = 1.726 ± 0.018 at z = 1091.3 (Komatsu et. al. 2010). 
Chi square in this case is denned as 

, _ (fl- 1.726) 2 
Xcmb - (aol8)2 (54) 

5.4 Joint Analysis with (H-z)+BAO+CMB 

Total chi-square function for our joint analysis: 

Xtot = Xohd + Xbao + Xcmb (55) 

The statistical analysis with Xtot gives the bounds on the 
model parameter specially on B. 



Model 




B 


A s 


a = 0.999,A.2V„ 


= 0.01 


0.00374357 


0.0593243 


a = 0.999,AiV„ 


= 2.00 


0.00434955 


0.0841938 


a = 0.001,AAT„ 


= 0.01 


0.00504535 


0.0400150 


a = 0.001,AJV„ 


= 2.00 


0.00600186 


0.0555941 



Table 2. Best-fit values for MCG: K = 1 

The contours between B and As for closed universe 
for q = 0.999 and a — 0.001 are shown in figs. (1) and 
(2) respectively. Figures 1(a), 2(a) are drawn for A_/V„=0.01 
and Figures 1(b) and 2(b) are drawn for A7V„=2.0. From 
fig. 1(a) which is plotted for a = 0.999 and AiV„ = .01 for 
closed universe it appears that the value of B lies in the 
range -0.05138 < B < 0.03351, -0.08234 < B < 0.08444, 
-0.1073 < B < 0.1274 at 68.3%, 95.4%, 99.73% confidence 
levels respectively. From fig. 1(b) which is plotted for a = 
0.999 and AiV„=2.0 for closed universe, the value of B lies 
in the range -0.05338 < B < 0.03251, -0.08434 < B < 
0.07745, -0.1103 < B < 0.1224 at 68.3%, 95.4%, 99.73% 
confidence levels respectively. It is evident from the analysis 



© 2012 RAS, MNRAS 000, 1-?? 



8 B. C. Paul and P. Thakur 



N 



a = .999 
AN^, .01 




a = .999 
AN^ 2.0 



>'ht 




[b] 

Figure 3. Constraints for open universe for a = 0.999 from 
OHD+SDSS+CMB Shift data 68.3%(Solid) 95.4% (Dashed) and 
99.73% (Dotted) contours. 
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Figure 4. Constraints for open universe for a = 0.001 from 
OHD+SDSS+CMB Shift data 68.3%(Solid) 95.4% (Dashed) and 
99.73% (Dotted) contours. 



that the range of values of B decreases with an increase of 
the effective neutrino parameter AN„. 

From fig. 2(a) which is plotted for a = 0.001 and 
ANu= 0.01 for closed universe, the value of B lies in the 
range -0.05042 < B < 0.03751, -0.09878 < B < 0.09906, 
-0.1274 < B < 0.1562 at 68.3%, 95.4%, 99.73% confidence 
levels respectively. From fig. 2(b) which is plotted for a = 
0.001 and AJV 1/ =2.0 for closed universe, the value of B lies in 
the range -0.05702 < B < 0.0419, -0.09988 < B < 0.1002, 
-0.1296 < B < 0.1595 at 68.3%, 95.4%, 99.73% confidence 
levels respectively. It is clear from the above analysis that 
the range of values of B increases with an increase in the 
effective neutrino parameter AN„. Thus for a given vlaue 
of a, the increase in neutrino parameter AN V increases the 
range of B parameter in the positive side. 

The contours between B and As for open universe for 
a = 0.999 and a = 0.001 are drawn in figs. (3) and (4) 
respectively. Figures 3 (a) and 4 (a) are drawn for AN„ 
=0.01 and Figures 3 (b) and 4(b) are drawn for AiV„=2.0. 
It is evident from fig. 3(a) which is plotted for a = 0.999 
and AN V = 0.01 for open universe, that the parameter B 
satisfies the following inequalities —0.05141 < B < 0.0178, 



Model 


B 




AJV„ 


a = 0.999 


0.007498 


0.107866 


0.100055 


a = 0.500 


0.010499 


0.110576 


0.100002 


a = 0.001 


0.016478 


0.114298 


0.100005 



Table 3. Best-fit values for MCG : K = -1 



Model B A s 



a 


= 0.999,AiV„ 


= 0.01 


0.00398596 


0.0629336 


a 


= 0.999,AjV„ 


= 2.00 


0.00782675 


0.133407 


a 


= 0.001,A7v„ 


= 0.01 


0.00528027 


0.0418902 


a 


= 0.001,AiV„ 


= 2.00 


0.00937448 


0.0817538 



Table 4. Best-fit values for MCG: K = -1 
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-0.08437 < B < 0.07932, -0.1063 < B < 0.1266 at 68.3%, 
95.4%, 99.73% confidence level respectively. In the fig. 3 (b), 
the allowed range of values of the parameter B for a = 
0.999 and AA jr „=2.0 for open universe are obtained which 
are given by -0.04921 < B < 0.02879, -0.08437 < B < 
0.07713, -0.1074 < B < 0.120 at 68.3%, 95.4%, 99.73% 
confidence levels respectively. It is evident that the domain 
of B decreases with an increase of the effective neutrino 
parameter AN„. 

Figure 4(a) is plotted for a = 0.001 and AiV„=.01 for 
open universe. In this case B lies in the following ranges 
-0.05262 < B < 0.03971, -0.09768 < B < 0.09686, 
-0.1274 < B < 0.1496 at 68.3% , 95.4% ,99.73% confi- 
dence levels respectively. From fig. 4(b) which is plotted for 
q = 0.001 and AiV„=2.0 for open universe, one can obtain 
viable cosmologies where B satisfies the following inequal- 
ities -0.05812 < B < 0.03751, -0.09768 < B < 0.09576, 
-0.1274 < B < 0.1584 at 68.3% , 95.4% , 99.73% confidence 
levels respectively. It is evident from the contours drawn in 
fig. (4) that the range of B increases with an increase in the 
effective neutrino parameter AN„ . But the positive range of 
values of B decreases. 



6 NUMERICAL ANALYSIS TO DETERMINE 
CONSTRAINTS ON THE EOS 
PARAMETERS IN BEYOND-DETAILED 
BALANCE SCENARIO 

In this section we use data to constrain the parameters of 
the MCG that were used in detailed balance scenario of the 
previous section. Stern data set for (H — z) data has been 
used along with BAO peak parameter and CMB shift pa- 
rameter. Chi-square minimisation technique has been used 
to determine the limiting values of the EOS parameters in 
the next subsections. 



6.1 (H-z) data as a tool for constraining 

The best-fitted parameters of the model considered here can 
be obtained by minimising the entity chi-square which is 
defined as 



X 2 ohd(H , Q Ko , As, B, a, f3, AN V , z) = (56) 
(H(H , n Ko ,A s , B, a, (3, AN v ,z) - H obs (z)) 2 



(57) 



where H obs (z) is the observed Hubble parameter at redshift 
(z) and a 2 z is the associated error with that particular ob- 
servation. Hubble parameter is given by 

(58) 



H(z) = H E(z) 
where we denote 

E 2 (z) = ^ibo(l+z) 3 +n co F(z)+n ro (l+z) i +n Ko (l+z) 2 

+[Qi + tt 3 (l + z) 4 + n 4 (l + z) 6 ] (59) 

with 

1-As 



F(z) = 



A s + 



a 3(l+B)(l+a) 



In this case H(z) — z data is taken from Stern Data analysis 
(Stern et. al. 2010). There are 12 data points of H(z) at 
redshift z used to constrain the MCG model. 



6.2 BAO peak parameter as a tool for 
constraining 

A model independent BAO (Baryon Acoustic Oscillation) 
peak parameter can be defined for low redshift (zi) mea- 
surements as: 



A = 



Jo 



E(z) 



Zl 



(61) 



where O m is the matter density parameter for the Universe. 
For a detailed description of the above defined parameter 
and related approximations reader is referred to ( Eisenstein 
et. al. 2005). The chi-square function can be defined as usual: 



2 

Xbao 



{A- 0.469) 2 
(0.017) 2 



(62) 



where we have used the measured value for A (0.469±. 0.017) 
as obtained by ( Eisenstein et. al. 2005) from the SDSS data 
for LRG (Luminous Red Galaxies) survey. 



6.3 CMB Shift Parameter as a tool for 
constraining 

Here the CMB shift parameter is defined as 



dz 

W) 



(63) 



where z ts is the z at last scattering. The WMAP7 data gives 
us R = 1.726 ± 0.018 at z = 1091.3 (Komatsu et. al. 2010). 
Chi square is defined as 



2 _ (R 

XCMB - 



1.726) 2 



(0.018) 2 



(64) 



(60) 



6.4 Joint Analysis with (H-z) +BAO+CMB 

We define total chi-square function for our joint analysis as: 
Xtot = Xohd + Xbao + Xcmb- (65) 

The statistical analysis with Xtot gives the bounds on the 
model parameter specially on B. 

Fig. 5(a) is plotted for a = 0.999 with best-fitted val- 
ues of f3, AN V and Qko- The parameter B satisfies the fol- 
lowing inequalities -0.05498 < B < 0.03808, -0.07989 < 
B < 0.07346, -0.1035 < B < 0.1128 at 68.3%, 95.4%, 
99.73% confidence levels respectively. Fig. 5(b) is plotted 
for a = 0.500 for best-fitted values of /3, AN„ and Q.ko- 
The parameter B in this case satisfies the folowing inequal- 
ities: -0.05741 < B < 0.04421, -0.08349 < B < 0.08287, 
-0.1060 < B < 0.1305 at 68.3% , 95.4%, 99.73% confidence 
levels respectively. Fig. 5(c) is plotted for a = 0.001 for best- 
fitted value of (3, AN V and Q,ko- We note that the parameter 
B satisfies the following inequalities -0.09257 < B < 0.0707 
, -0.1326 < B < 0.1493, -0.1727 < B < 0.2247 at 68.3%, 
95.4%, 99.73% confidence levels respectively. It is evident 
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Figure 5. Constraints in beyond detailed balance for a = 
0.999, a = 0.500,a = 0.001,from OHD+SDSS+CMB Shift data 
68.3%(Solid) 95.4% (Dashed) and 99.73% (Dotted) contours. 



[b] 
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Figure 6. Constraints in beyond-detailed balance for a = 0.999 
from OHD+SDSS+CMB Shift data 68.3%(Solid) 95.4% (Dashed) 
and 99.73% (Dotted) contours. 



that the allowed range of values of the parameter B be- 
comes larger compared to that of the detailed balance sce- 
nario (Paul et. al. 2012). 

Fig. 6(a) is plotted for a = 0.999 and A7V„=0.01 for 
best-fitted value of j3 and Qko, it is evident that B can 
take any value in the following ranges —0.05338 < B < 
0.03351, -0.07935 < B < 0.07445, -0.1013 < B < 0.1224 
at 68.3%, 95.4%, 99.73% confidence levels respectively. Fig. 
6(b) is plotted for a = 0.999 and AA^^.O, it is evident that 
the value of B lies in the range -0.05462 < B < 0.03617, 
-0.07978 < B < 0.07662, -0.1032 < B < 0.1162 at 68.3%, 
95.4%, 99.73% confidence levels respectively. The figs. 6(a)- 
6(b) show that the range of permissible values of B decreases 
with an increase in the effective neutrino parameter. 

Fig. 7(a) is plotted for a = 0.001 and A7V„=0.01 for 
best-fitted value of f3 and Q.ko, it is evident that the per- 
missible values of B now lies in the range —0.06682 < B < 
0.0407, -0.09434 < B < 0.08918, -0.1206 < B < 0.1416 
at 68.3%, 95.4%, 99.73% confidence levels respectively. Fig. 
7(b) is plotted for a = 0.001 and A7V„=2.0, it is evident that 
the values of B lies in the range -0.06347 < B < 0.04844, 
-0.09244 < B < 0.09241, -0.1194 < B < 0.1414 at 68.3%, 
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Figure 8. Equation of state parameter in beyond-detailed bal- 
ance scenario 



1 2 

Ptot = Pc + ~pr H X 



K 

Aa 4 



3A 



Ptot = pc + Pb + Pr H J 



3K 2 
Aa 4 



+ 3A 



(66) 



(67) 



In the case of beyond-detailed balance scenario the total 
pressure and the energy density is given respectively as 



[b] 

Figure 7. Constraints in beyond-detailed balance for a = 0.001 
from OHD+SDSS+CMB Shift data 68.3%(Solid) 95.4% (Dashed) 
and 99.73% (Dotted) contours. 
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Pc + gPr 
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(68) 
(69) 



95.4%, 99.73% confidence levels respectively. The contours 
drawn in figs 7(a) and 7(b) show that the range of permissi- 
ble values of B now decreases with an increase in the effec- 
tive neutrino parameter. We note that the allowed range of 
values of the parameter B, decreased appreciably here com- 
pared to that obtained from figs. 5 (a)-5 (c). This signifies 
the fact that as the contribution of dark radiation increases 
(through effective neutrino parameter) the range of admis- 
sible values of B decreases in the case of beyond detailed- 
balance scenario which is same as one obtains in the case of 
detailed balance scenario. 



7 VIABILITY OF MCG IN HL GRAVITY 

In this section we discuss some of the implications of the 
present scenario. Here we determine the evolution of the 
equation of state parameter of the total cosmic fluid of the 
universe which is defined as w(z) — SML with the total 
pressure and energy density in the case of detailed-balance 
scenario. The pressure and energy density are given by 



Here we replace the scale factor by redshift parameter, and 
the expression for density parameter and the Hubble param- 
eter are expressed in terms of the redshift parameter in the 
equation of state. Consequently we get 

„,( z ) = 2*2*. (70) 

ptot 

From the plot of w(z) with z, for beyond detailed- 
balance scenario we note that at high redshift (i.e., early 
times) it attains a fixed value | since radiation dominates 
in that epoch. In the intermediate redshift it behaves as dust 
for quite a long time. It is observed that the equation of state 
parameter picks up negative values at small redshift, i.e., at 
very recent past. In the case of closed or open universe the 
present value of the equastion of state parameter is found to 
be negative (-0.7) which admits a late accelerating universe. 

In order to check the validity of the scenario we em- 
ployed the best-fit values of the parameters of the MCG 
to find supernovae magnitudes (/i) at different redshift (z) 
and plotted n vs. z curve. We compared these with original 
curves of Union2 data and observed an excellent agreement. 
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Figure 9. The comparison of the Union2 data with the best-fit 
values in beyond-detailed balance 



8 DISCUSSION 

In this paper we present cosmologies with modified Chap- 
lygin gas (MCG) in HL gravity scenario taking into ac- 
count detailed balance and beyond detailed balance condi- 
tions both in the presence and absence of dark radiation. 
The equation of state of MCG has three unknown param- 
eters. The permissible values of these parameters are ex- 
plored from the observed data. Using data from different 
observations, namely, H(z) — z (OHD), BAO peak, CMB 
shift parameter data, we determine the admissible values of 
the EOS parameters considered here. In the MCG we have 
a parameter B that represents the matter part. We analyze 
and determine the allowed range of values of B for viable 
cosmologies. The analysis is carried out here both in the 
case of open and closed universe at 68.3%, 95.4%, 99.73% 
confidence levels. 

In a close universe we note that the permissible values 
of B parameter lies in the range -0.05702 < B < 0.0419, 
-0.09988 < B < 0.1002, -0.1296 < B < 0.1595 at 68.3%, 
95.4%, 99.73% confidence levels respectively for a maximum 
value of effective neutrino parameter. The range of B ob- 
tained here is less than that obtained for best-fit value of ef- 
fective neutrino parameter when dark radiation is not taken 
into account (Paul et. al. 2012). 

It is evident from the figures 2 (a) and 2(b) that the 
range of B increases with an increase in the effective neutrino 
parameter A7V„. In an open universe, the parameter B lies 
in the range -0.05812 < B < 0.03751, -0.09768 < B < 
0.09576, -0.1274 < B < 0.1584 at 68.3%, 95.4%, 99.73% 
confidence levels respectively for maximum value of effective 
neutrino parameter. The domain of B is found to be less 
than that obtained for best-fit value of effective neutrino 
parameter without dark radiation (Paul et. al. 2012). 

It appears from the analysis of a close and an open uni- 
verse that in a close universe the domain of admissible values 
B is comparatively narrower than that of an open universe 
at 68.3%, 95.4%, 99.73% confidence levels respectively. We 
alaso note that the B may take negative values in this case. 
The negative value of B implies existence of exotic matter. 
In the literatures (Fabris et. al. 2010) the acceptable value 



of B as was predicted to be very small, which once again 
gets support from our analysis. 

In the beyond-detailed balance scenario there are six 
free parameters, namely Qko , As, B, a, f3, AN V . It is found 
that the entire range of effective neutrino parameter is con- 
sistent with observations from our numerical analysis. The 
dependence of the extreme values of the neutrino parameter 
on other parameters are also shown in figs. (6) and (7). 

The contours drawn in fig. (5) for beyond-detailed bal- 
ance with different a for best-fitted value of f3, AiV„ and 
Qko projects the admissible values of B which lies in the 
range -0.09257 < B < 0.0707, -0.1326 < B < 0.1493, 
-0.1727 < B < 0.2247 at 68.3%, 95.4%, 99.73% confidence 
levels respectively. Thus the range of B in this case become 
larger than that of detailed balance scenario without dark 
radiation (Paul et. al. 2012). 

The contours drawn in fig. 6(a) which is plotted for 
a = 0.999 and A7V„=0.01 for best-fitted value of j3 and 
f^Ko, projects the admissible values of B which lies in the 
range -0.05338 < B < 0.03351, -0.07935 < B < 0.07445, 
-0.1013 < B < 0.1224 at 68.3%, 95.4%, 99.73% confidence 
levels respectively. The contours drawn in fig. 6(b) which is 
plotted for a — 0.999 and AiV„=2.0, permits the values of 
the parameter B which lies in the range —0.05462 < B < 
0.03617, -0.07978 < B < 0.07662, -0.1032 < B < 0.1162 
at 68.3%, 95.4%, 99.73% confidence levels respectively. It is 
clear that the range of values of B decreases with an increase 
of the effective neutrino parameter. 

The contours drawn in fig. 7(a) which is plotted for 
a = 0.001 and AiV„=0.01 for best-fitted value of (3 and 
Qko gives the allowed values of B which lies in the range 
-0.06682 < B < 0.0407, -0.09434 < B < 0.08918, 
-0.1206 < B < 0.1416 at 68.3%, 95.4%, 99.73% confidence 
levels respectively. The contours drawn in fig. 7(b) which is 
plotted for a = 0.001 and AiV„=2.0 gives the allowed val- 
ues of B which lies in the range -0.06347 < B < 0.04844, 
-0.09244 < B < 0.09241, -0.1194 < B < 0.1414 at 68.3%, 
95.4%, 99.73% confidence levels respectively. It is clearly vis- 
ible from above analysis that the range of values of B de- 
creases with an increase in effective neutrino parameter. 

We note that the range of values of B decreases appre- 
ciably here compared to that obtained from the figs. 5 (a)- 5 
(c). This signifies that as the contribution of dark radiation 
increases (through effective neutrino parameter) the contri- 
bution to the permissible range of values of B decreases in 
beyond detailed-balance scenario like that one obtains in the 
case of detailed balance scenario (Paul et. al. 2012). In figure 
(8) we plot the variation of the total equation of state pa- 
rameter w(z) with redshift parameter z for beyond detailed- 
balance scenario. The curve shows the evolutionary phases 
of the universe efficiently. It is evident that at high redshift 
(early times) the of equation of state parameter attains j , in- 
dicating radiation domination in that epoch. However, in the 
intermediate redshift we note that dust dominates through 
MCG for quite a long period of time. The best fit values of 
B, As and AiV„ are shown in tables (1) and (2) for close 
universe and tables (3) and (4) are shown for open universe. 

Using the best-fit values in beyond-detailed balance sce- 
nario jU vs redshift curve is plotted in fig. (9) and the figure 
is compared with Union Compilation data. It is evident from 
the figure that cosmologies in Horava-Lifshitz gravity with 
MCG fits well with the experimental result. 
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In this analysis we studied the dependance of extreme 
values of AN V on other parameters in detailed balance sce- 
nario with MCG. We also present here results obtained on 
the beyond-detailed balance scenario. In the case of beyond- 
detailed balance scenario (BDB) there are two more free 
parameters and found that the theory is rich and all the 
features of cosmologies can be accommodated with MCG. 
Earlier in the Einstein-frame, MCG is employed to obtain 
viable cosmological models (Thakur, Chose & Paul 2009; 
Jianbo et al. 2010). Here MCG is employed in the HL grav- 
ity and determined various physical parameters of the uni- 
verse which gets supports from observations. However, the 
present analysis does not enlighten the conceptual issues in 
HL gravity. It is important to look into details why the neu- 
trino parameter is very small in HL gravity with MCG which 
will be discussed elsewhere. 
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